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Abstract. We study the problem of charging a quantum battery in finite time. We 
demonstrate an analytical optimal protocol for the case of a single qubit. Extending 
this analysis to an array of N qubits, we demonstrate that an N-fold advantage in 
power per qubit can be achieved when global operations are permitted. The exemplary 
analytic argument for this quantum advantage in the charging power is backed up 
by numerical analysis using optimal control techniques. It is demonstrated that the 
quantum advantage for power holds when, with cyclic operation in mind, initial and 
final states are required to be separable. 
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1. Introduction 

The maximum work that can be extracted unitarily from a given quantum state p with 
respect to a reference Hamiltonian Hq was called ergotropy by Allahverdyan, Balian, 
and Nieuwenhuizen (Allahverdyan et ah 2004). In this extraction scheme the final 
system state is unique up to degeneracies at the level of the Hamiltonian and belongs 
to the class of passive states (Pusz & Woronowicz 1978, Lenard 1978). Such states are 
characterised by being diagonal in the energy eigenbasis of the Hamiltonian and ordered 
with decreasing eigenvalues corresponding to increasing energies. 

This paradigm of unitary work extraction has recently been extended to scenarios 
where multi-partite systems are considered. It has been shown by Alicki and Fannes 
(Alicki & Fannes 2013) that the extractable work can be increased by allowing for 
entangling operations. Subsequently, Hovhannisyan et ah (Hovhannisyan et al. 2013) 
demonstrated that, irrespective of an entangling operation being required for optimal 
work extraction, the entanglement created in the extraction processes can indeed be 
minimal and even vanish. They argue that the cost of this minimisation is a longer 
process duration and consequently conjecture that there may be a relation between the 
power of such a process and the entanglement created. This line of reasoning has been 
extended to explore the work cost in the creation or erasure of more general correlations 
(Giorgi & Campbell 2015, Huber et al. 2015, Perarnau-Llobet et al. 2014, Bruschi 
et al. 2015). In this article we are interested in the inverse process of work extraction: 
charging a quantum battery. 

We are motivated by the following question: Are there truly quantum effects 
or phenomena intrinsic to a specific thermodynamic process that offer an operational 
advantage over their classical counterpart? After all thermodynamics should be 
sensitive to the underlying microscopic description. For example, the efficiency of 
an engine should always be limited by the Carnot bound irrespective of whether the 
working medium is comprised of quantum or classical components. However, for non- 
equilibrium protocols the question is more subtle. For instance very recently it has been 
demonstrated that it is possible to use non-equilibrium short cuts to adiabaticity to 
boost the power of engine cycles without compromising efficiency (Deng et al. 2013, del 
Campo et al. 2014). One may then wonder whether it is in the finite time operation of 
devices that quantum mechanics offers an advantage. The present article provides an 
intuitive example of one such process where quantum correlations provide an advantage: 
When taking the finite time nature of a charging process into account. In this case we 
demonstrate that the non trivial quantum nature of both the Hamiltonian and the state 
leads to a substantial advantage. That is, entangling processes, which allow for the 
creation of quantum correlations, can achieve a higher power than local operations. 
This is the case even when correlations are only allowed to appear during the process 
but not in the input and output states. 

This paper is structured as follows. In Sec. we outline the general problem of work 
extraction and powerful driving in a quantum system. We then derive a power-optimal 
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driving for a single qubit - a case where constraints on the driving are unambiguous - 
in Sec. Sec. extends the results to an array of N qubits and we demonstrate an 
advantage in power that stems from the permittance of entangling operations both for 
an initial pure and an initial thermal state. We conclude in Sec. 

2. Charging a quantum battery 

A battery is a physical system that stores energy. In this paper we use quantum systems, 
which store energy in the energy levels and coherences, for this purpose (Alicki & 
Fannes 2013). The internal energy of a quantum system (the battery in this case) is then 
given by tr[p.^o]) where p is the state of the battery and IIq is its internal Hamiltonian. 
Charging a battery is to change its state from p to a more energetic state p' such that 
tr[(p' — p)Ho] > 0. Conversely, using the battery will take it a lower energy state p" 
such that tr[(/5" — p)Hq] < 0. We restrict both charging and discharging processes to 
be cyclic unitary processes. This is achieved by applying an external time-dependent 
potential V{t) for time T, in addition to the internal Hamiltonian of the battery Hq. As 
a consequence of unitarity the spectrum {pi} of any accessible battery state is fixed [f} 
The lowest energy state, vr, is called a passive state, and analogously the highest 
energy state, a), a maximally active state; both of them are defined with respect to Hq. 
Let us express the internal Hamiltonian with increasing energy levels as 

Ho := \ei) (ej| with < e^+i. (1) 

The passive and active states are respectively 

^ Pi |c) (cl with Pi > Pi+i and (2) 

i 

to :='^Pi\ei) {€i\ with Pi < Pi+i. (3) 

i 

Note that all thermal states are passive. 

If the battery is in a generic state p = '^iPi |rj) {rf (with Pi > Pi+i), we can use 
all available energy and the battery will end up in its passive state. The amount of 
extractable energy from this battery is called ergotropy >V. It is defined as the cyclic 
work extractable with respect to the internal Hamiltonian, W := tv{Hop) — tr(iLo^)- 
Once a quantum state p has been transformed to a passive state tt, no more work 
is extractable from that state via cyclic unitary transformations. Therefore ergotropy 
quantifies the available energy in a battery. Conversely, a fully charged battery is in 
state uj and the system cannot be charged further via a cyclic unitary process. Note 
that states tt, u, and p are unitarily related to each other. 

In practice it is often desirable to charge a battery quickly, or better: with maximal 
power. Let us consequently address the question of powerful charging of quantum 

I When comparing two different spectra {qi} and {pi}, we note that {pi} allows for both a lower 
minimum-energy state as also a higher maximum-energy state if it majorises {qt} (Marshall et al. 2011). 
Consequently, batteries in a pure state have higher capacity than ones in mixed states. 
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batteries. For the charging protocol we start with a generic state p and apply cyclic, 
unitary evolution for some optimal time T, 

Ht = Hq + Vt with td = 0, for t < 0 and t > T, (4) 

where the index t indicates time-dependence. The protocol duration T is not fixed but 
rather part of the optimal solution. Here we are not particularly interested in the specific 
output state p' but rather in the average work fW) done in the process and, importantly, 
the average power (P) = (W)/T. Alternatively, one could study other process-duration 
dependent quantities. For instance, maximising (P) x (kF) will guarantee that power is 
not optimised at the cost of average work. In order to be fully general we thus consider 
the family of objective functions 

T {Pr {W)'-"-= , (5) 

where T is the process duration and 0 < a < 1 remains a free parameter. For the 
extremal case of a = 1 the expression for power is recovered. Optimisation for a = 0, 
on the other hand, recovers the expression for ergotropy as derived in (Allahverdyan 
et al. 2004) - This will be demonstrated in a forthcoming publication (Vinjanampathy 
et al. 2015) as part of a general optimal control approach to extremising P. 

A few remarks are in order: 

(1) Cyclicity of the process is an important requirement since we are interested in 
the power. Consider a non-cyclic process such as a sudden quench, Hq —> H', which 
does a finite amount of work on the system in an instant and therefore produces infinite 
power. Such a scenario renders the problem of studying power trivial. If, on the other 
hand, the process is required to be cyclic then we must have Hq ^ H' ^ Hq, and, for a 
sudden quench, the objective function vanishes both in the numerator and denominator. 

(2) One could imagine driving the system infinitely fast between two given states 
and let the denominator of Eq. approach zero. Such infinitely powerful driving, 
however, would be in conflict with so-called quantum speed limits which bound the 
evolution time by the inverse mean energy (Margolus & Levitin 1998), or, alternatively, 
its variance (Mandelstam & Tamm 1945) and apply to time-dependent unitary evolution 
between any two quantum states (Deffner & Lutz 2013). These speed limits are not 
prescriptive for determining optimal driving because they depend on the the system 
state at each instant of time. Alternative bounds, however, have been developed that 
depend on norms of the driving Hamiltonian (Uzdin et al. 2013). 

(3) Such bounds on the Hamiltonian are natural constraints that must be 
incorporated into the analysis of quantum batteries operating at finite time. In the 
following sections, we will propose bounds that are external constraints in the sense that 
they depend on the driving Hamiltonian alone but not the system state. In particular, 
we will restrict the trace norm of the driving Hamiltonian ||.Hf|| < Emax with the gauge 
convention that its lowest eigenvalue is zero. This corresponds to restricting the maximal 
energy available for external driving. 
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Let us now derive an optimal driving with respect to the objective function, Eq. 
for a single qubit. Without loss of generality we take Hq = |1} (1|. We parametrise 
the driving Hamiltonian Ht with control functions vf,v^,v^ for the Pauli operators 
ax, ay, dz- The three control functions can be grouped in a vector Vj = (uf, uf, uf and 
similarly for the Pauli operators <t = (da,, ay, az)^. The instantaneous eigenvalues of Vt 

are = ±|lit| = A thermodynamically sensible constraint is achieved 

by bounding the difference between the instantaneous eigenvalues by a maximum value 
Emax ~ that is, Xf — Xf < Emax for all times t in accordance with the trace norm 
introduced above. 

We now follow recent results, where full analysis of a general time-optimal control of 
qubit operations was performed (Hegerfeldt 2014, Hegerfeldt 2013). More generally this 
is related to the so-called quantum brachistochrone problem (Carlini et al. 2003). These 
approaches differ from the present analysis in the sense that they are concerned with a 
process which takes some input state p to a specific state p' whereas we are concerned 
with the average work that is done in a driving process. Whilst we are interested in 
preserving the (qubit) state’s purity it makes no difference what phase it has. In the 
Bloch Sphere picture, we are only interested in the state’s 2 : coordinate. In this context, 
p can be any qubit state (active or passive), i.e., parametrised by angles 9 and f, and 
a radius r which remains constant during evolution. The value of f has no bearing on 
the result. As will be seen below this leads to zero driving along in the optimal case. 

In order to optimise the power we now invoke the von Neumann equation for the 
state’s unitary evolution: 

= [Ho + Vt, Pt] = 1 + at.d-], (6) 

where pt is parametrised by its Cartesian decomposition in terms of Pauli operators, 
i.e. a.t = {af,a!j.,alY' = r{sm6tCos4>t, sindfsin^j, cos 9t)^. With Einstein summation 
convention and using the Levi-Civita symbol the evolution becomes 
d 


dt 


pt = vla^ejkidi. 


(7) 


We want to achieve maximum average power (P) = (IT) /T (or indeed T = (IT) /T") 
over the duration T of the whole process. To reflect cyclicity, work (and consequently: 
power) is here defined with respect to the initial Hamiltonian Hq: Charging increases 
the energy in the state pt with respect to the time-independent reference Hq. Optimality 
is hence achieved by first optimising E^^r[ptHo] at each instant in time over the driving 
Vt and then finding the optimal process duration T. Using [Equation 7] we obtain 


tr 




= [vl cos (ft — vf sin (f)t)r sin 9t. 

The optimal protocol under the above constraint is found for 


— En 


:Sin ft. 


— E 

^max 


cos ft, 


Et = 0 . 


( 8 ) 

(9) 

( 10 ) 
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This implies that in the optimal case no driving happens in the direction of the reference 
Hamiltonian. The solution corresponds, unsurprisingly perhaps, to driving along the 
geodesic with hxed r and ft = fo at constant angular speed E^ax, that is 

Ot = Oq + Emaxt- ( 11 ) 

For a fixed input state we now want to optimise the average power, or more generally 
the function 


(W) ^ tr|gop'| - tr[gop] 

'Y'a 'Y'a 

f Emax V / cos 6*0 COS 6 *t \ 
\6*t — 6*0/ \ 2 2 / 

[cos 6*0 - cos( 6 *o + EmaxT)]. 


( 12 ) 

(13) 

(14) 


Linearity of the driving time in Ot (cf. [Equation 11 ) allows us to choose either variable 
for optimisation. The optimal solution is found for the driving time Tm that satisfies: 


cos(0o T E^g^fPjn) 
This equation can be 


E T 


sin( 


(^0 T E^axTm) COS^o- 
a 

understood geometrically by rewriting it as 


(15) 


(16) 

a 

which is fulfilled when the qubit has travelled a vertical distance /S.zt on the Bloch 
Sphere which equals the length of its projection onto the a:-?/-plane, , weighed by a 
factor EmaxT/a. 

The solution for T has a non-trivial maximum which is different for different a 
as can be seen for the example in Fig. |l(a)[ The cone of optimal output states for 
initially thermal states is pictured on the Bloch Sphere in Fig. |l(b)[ We point out that 
for certain initial conditions it is possible to achieve finite power, or indeed W, at time 
T —>■ O"*". This can be regarded as pathological and it is worth pointing out that for 
a < 1 and 9 < % the optimum is always reached for finite T. 


4. Role of entanglement in charging an array of quantum batteries 


Having found the maximally powerful evolution of a single qubit we now proceed by 
looking at larger systems, in particular arrays of N qubits, i.e., as shown in 

Fig.i The aim of the following section is to demonstrate that when global, entangling 
operations are allowed on the array (rather than local, parallel operations for each qubit) 
an increase in power linear in N can be achieved. 


§ For a given input state Equation 11 already lets us answer two related questions: What is the 
maximum work that can be achieved in time T? What is the minimum time necessary to achieve work 
W (parametrised by 0)7 
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Figure 1. (Colour online.) (a) Taking the generic, thermal input state f = 
|[(1 ± r) |0) (0| + (1 =F ?’) |1) (1|] average power (P) := [a = 1, red, solid], and 
[a = |, blue, dashed] are plotted against T. Maximal power is reached for 
9f « 0.747r in this case, (P) (W) is maximal at 9f « O.SQtt [grey, dotted lines]. The plot 
was created with Emax = 1 for simplicity and the functional values are given in units 
of 1/r. (b) For a battery in an initially thermal state (i.e. 9q = tt, thick black line), 
and Emax = this graph shows the cone of states that achieve maximum average 
power on the Bloch Sphere. Maximal power is reached for 9t ~ 0.747r. Note that a 
cone only appears for initial states on the 2 ;-axis. For all other states (p is fixed and 
the optimal output state lie on a ray at <f>T = 4>o- The green line depicts an exemplary 
trajectory for r = | and ^ = 0. 


The reference Hamiltonian is the sum of the local Hamiltonians from the single qubit 
case = Ylik |l)fc (l|fc where k labels the kth qubit. For now, we don’t want 

to decrease the purity of the local states (with cyclic operation in mind). We define this 
- the condition that there is no reduction in the purity of the marginals at the end of 
the protocol - as the condition of no degradation. Restricting the analysis to pure states 
for now, we are interested in going from := |0)®'^ to (a |0) + 611))®-^. Since a 

full optimisation over all possible unitary transformations is not analytically feasible we 
will contend ourselves with finding a specific evolution that shows an improvement in 
power per qubit compared to the single qubit case. 

Let us consider full charging such that := For parallel 

driving, e.g. Eq. HI the time-independent driving Hamiltonian during 0 < t < T is 

N N 

Emaxicos - siu (g) ■ (17) 

k j^k 

It possesses A/^ + 1 equidistant eigenvalues with the gap between the largest and the 
lowest eigenvalue given by NEmax- In order to allow for a fair comparison we will hence 
equally allow for a maximum energy gap of Emax = NEmax in the case of global driving. 
The work done is trivially (IF) = N - that is, the average work per qubit is unity and 
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Figure 2. (colour online) Consider an array of qubits (A) that can either be driven 
in parallel (B), for instance as prescribed by Eq. 17, or globally (C), i.e. as given by 
Eq.lTsl 


the driving time is Tpar = T^/Emax- Global driving is achieved by evolution with[j|] 

= BSd"’ - EZl dl'"’) (0^1 + lo'"') <1'~'|) (18) 

for a duration T = The average work in this process is again (W) = N. The 

average power consequently reads 

Z7'(^) 

-^max 

TT 

and depends on the driving constraint EmJx for . Permitting for the same constraint 
as in the parallel case, that is Emax = NE^ax, an Af-fold increase in power per qubit 
is achieved. Interestingly, the advantage from correlations appears even though initial 
and final states are separable. 

The quantum speedup observed here can be understood in part as a consequence 
of the shorter distance that has to be travelled through state space when entangling 
operations are permitted. To calculate and compare the path length in state space 
between the global and local pure state case (i.e., —?• we first note that in 

the global case the proposed evolution does in fact prescribe a path along a geodesic. 
As a consequence the path length is directly given by the Bures angle C between the 
input and the output state. For pure states the Bures angle is equal to the Fubini-Study 
distance and we directly see (Bengtsson & Zyczkowski 2006) 

£(|oOT).|im» 

= arccos |(0<«l|lW>|=^ (20) 

Li 

II This global driving is in fact optimal for the current scenario. Generally, global driving between any 
two states |a) and \b) is optimised by evolution with oc |a) {b\ + \b) {a\ where the proportionality is 
given by the external constraint. 
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In the case of parallel driving we note that the sqnare of the line element is additive 
under tensor product. That is, for Ifjglobal) = ds^giohai ~ Accordingly 

the path length in this case scales with '/N. 



Figure 3. (colour online) Numerical optimization results for minimum time t± for 
complete charging of N qubits, where the state of the battery rotates from 
to The maximum eigenvalue of the Hamiltonian is scales linearly with the 

number of qubits, and we note that the full charging time tj_ is commensurate with 
1/N. The solid line represents using the single qubit charging protocol for all N qubits. 
The dashed line represents the 1/N line, whereas the circles represent the solution to 
numerical optimisation. The entanglement dynamics for the four qubit numerical 
optimization is shown in (b), by plotting the entropy of the reduced state obtained 
by tracing over the last two qubits. We note that at the beginning and at the end of 
the dynamics, the qubits are disentangled, indicated by vanishing entropies. For the 
numerical optimisation the eigenvalues of the instantaneous Hamiltonian were bounded 
between 0 and Emdx = 1. Since the absolute scaling E^ax of the Hamiltonian is 
arbitrary, the single qubit time serves as a normalization to see the behaviour exhibited. 

For comparison, let us numerically study the same charging process, i.e., namely 
going from ^ Fig. I shows results of a numerical optimisation of full 

charging between the number of qubits in (a), plotted on the a:-axis, and the time 
that the quantum system takes to rotate between the states. The minimum time 
represents an optimisation over all possible time-independent Hamiltonians such that 
their eigenvalues are bounded by the same constraint as above (a factor linear in N). 
We see that the charging time is commensurate with 1/N behaviour. We note that 
the dots represent the numerical optimisation of full charging, whereas the dashed line 
represents the 1/N behaviour, as a guide to the eye. The solid line represents parallel 
charging of all N qubits using the single qubit protocol. Finally, the plot in (b) represents 
the evolution of the entanglement of the quantum system for the four qubit numerical 
solution. Plotted is the reduced entropy of the first two qubits, showing the generation 
of entanglement across the 12-34 bipartition. Note that as the system approaches the 
final state, corresponding to the optimal value of power, the qubits disentangle. This 
is commensurate with the demand that all batteries be available for individual use at 
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the final time with loss of purity. This exemplifies the role of entanglement during the 
protocol as being responsible for the quantum speedup. 

5. Conclusion 

In this article we have given a full derivation of a qubit protocol that achieves maximum 
power when charging a quantum battery - a work qubit - under constrained driving. 
With cyclic operation in mind we allowed for unitary evolution, hence keeping the state’s 
purity constant. In extension we then examined the charging of an array of N work 
qubits. Using a specific exemplary evolution we demonstrated Wfold advantage in 
power per work qubit. This example was presented under the operational constraint 
that the purity of the state must again be conserved. This highlights the potential for 
quantum enhancement of devices working under non equilibrium conditions. For future 
work it will be interesting to extend this scenario and allow for degradation - that is, 
an extra purifying stage is added after the discharging of the battery array to make 
up for any loss in purity during the global charging process. Furthermore, in order to 
allow for realistic noise and decoherence processes it will be of paramount importance 
to generalise the approach to open systems (Binder et al. 2015). We hope that this work 
will inspire further investigations into the advantages of quantum effects in finite-time 
thermodynamic processes. 
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